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Abstract 

This paper is devoted to a special class of the Bienayme-Galton- Watson pro- 
cesses with a countable type space whose probability generating functions 
are linear-fractional. For such processes using various tools (contour process, 
spinal representation, Perron-Frobenius theorem for countable matrices, re- 
newal theory) we thoroughly investigate the case of i?-positive recurrence 
with respect to the type space. 
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1. Introduction 

Branching processes is a steadily growing body of mathematical research 
having applications in various areas, primarily in theoretical population bi- 
ology [13], [16], [19]. A basic version of branching processes, called the Bi- 
enayme-Galton- Watson (BGW) process, describes populations of particles 
which live one unit of time and at the moment of death give birth to a 
random number of new particles independently of each other. In the sin- 
gle type setting the consecutive population sizes {Z^'^^n>o form a Markov 
chain with the state space {0,1,2,...}. An important analytical tool for 
studying branching processes is the probability generating functions. Given 
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(f)(s) = E(s^'^^), the n-th generation's size is charaterized by the n-fold iter- 
ation of (/){■) 

E(s^'"') = 0(...(0(.))...). 

Here and elsewhere in this paper we always assume that a branching process 
starts from a single particle. 

The case of hnear-fr actional generating functions is of special interest 
as their iterations are again linear-fractional functions allowing for explicit 
calculations of various entities of importance (see [3], p. 7). Such explicit 
results, although being specific, illuminate the known asymptotic results con- 
cerning more general branching processes, and, on the other hand, may bring 
insight to less investigated aspects of the theory of branching processes. 

In the multi-type setting particles still reproduce independently but now 
the number of offspring may depend on the mother's type. A flexible family of 
population models is obtained by means of BGW-processes with countably 
many types ([4], [H], [16], [23]). These are infinitely dimensional Markov 
chains 

Z(") = (Z;"),4"\...), n = 0,1,2,..., 

whose i-th component Z^"'^ gives the number of particles of type i existing 
at time n. In this paper we study the class of such branching processes with 
the generating functions for vectors Z*^") all being linear-fractional. As shown 
in Section |2] a linear-fractional BGW-process with countably many types is 
fully specified by a triplet of parameters (H,g,m), where H = {hij)fj^^ is 
a sub-stochastic matrix, g = [gi, g2, . . .) is a proper probability distribution, 
and m is a positive constant. 

A linear-fractional branching process with parameters (H, g, m) has the 
following reproduction law. A particle of type i has no offspring with proba- 
bility hio = 1 — hij. Given that this particle has at least one offspring, 
the type its first daughter is j with probability hij/ (1 — hio), and the number 
of subsequent daughters has a geometric distribution with mean m. With the 
exception of the first daughter the types of all other offspring particles are 
chosen independently of mother's type according to the common distribution 

g- 

Theorem |3] in Section [3] states that the probability generating function 
for the n-th generation of the branching process with parameters (H, g, m) 
is linear-fractional. This theorem presents a general formula for the gener- 
ating function of Z*^"). An important step in obtaining this formula uses a 
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spinal representation argument making the derivation more transparent as 
compared to the argument used in [12] for the finite-dimensional case. 

Another way to treat a BGW-process with countably many types, which 
is mentioned in Section |4], is to see the sequence of first born daughters of 
a particle as a single individual evolving in time by changing its type. We 
emphasize these two different interpretations of the same branching system 
by referring to the reproducing entities either as particles living one unit of 
time and replaced by random numbers of offspring (non-overlapping gener- 
ations) or individuals with random life lengths which are able to produce 
random numbers of offspring at different times along their life span (overlap- 
ping generations). In the linear-fractional case the evolution of individual's 
type is governed by the sub-stochastic matrix H. Every time unit during 
its life span, except the moments of its birth and death, the individual gives 
birth to an independent geometric number of offspring with mean m and the 
initial type for each offspring is chosen according to the distribution g. 

Assume that the type of the initial particle has distribution g. Then the 
total population size Z*^"^ = Z*^"-*!* of the linear-fractional BGW-process has 
the mean M„ := E(Z'-"^) which is computed as 

Mn = gM"l\ (1) 

where 

M = (m,,)^^^,, m,, = E(Z]')|Z(°) = e,) (2) 

is the matrix of the offspring means. The sequence {Z(")}„>o forms a single- 
type Crump-Mode-Jagers (CMJ) process, which we will call a linear-fractional 
CMJ-process. This CMJ-model is not restrictive about the individual's life 
length L > 1 distribution (see Example 8 in Section |4]) but individual's 
reproduction point process must follow a very specific pattern: at times 
1, . . . , L — 1 the individual produces independent and identically distributed 
geometric numbers of daughters. Making birth events in the linear-fractional 
CMJ-processes very rare (by choosing the row sums of H to be close to 
zero) and rescaling time accordingly we arrive at the family of continuous 
time CMJ-processes with a Poisson point process reproduction. This family 
of branching processes was investigated in [17j emphasizing the role of the 
so-called contour processes of the corresponding planar genealogical trees. 

In Section [5] we remind the concept of a contour process generated by a 
planar BGW-tree. In the multivariate linear-fractional framework a jumping 
version of the contour process (in the spirit of [17]) has a nice Markovian 
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structure of a constant speed descent with independent and identically dis- 
tributed upward jumps until absorbtion at the level —1. We apply the con- 
tour process method to prove that given the reproduction parameters (g, m) 
are independent on mother's type, the vectors Z^") have also linear-fractional 
distributions for all times n as stated in Theorem [3l 

Section |6] introduces a double classification of the branching processes 
based on the Perron-Frobenius theorem for countable matrices. Besides the 
usual classification into subcritical, critical, and supercritical branching pro- 
cesses in the case of infinitely many types one has to distinguish among 
i?-transient, i?-null recurrent, and i?-positively recurrent cases depending on 
the corresponding property of the mean matrix (j2]). 

The main result of this paper. Theorem 10 in Section [7| among other 
statements contains a transparent criterion for /^-positive recurrence of the 
linear-fractional branching processes. In the language of the CMJ-processes 
this criterion requires that the corresponding Malthusian parameter is well 



defined and the mean age at childbearing is finite. Theorem 10 is proven in 
Sections H] and [H 

In Section 10 we present the basic asymptotic results for subcritical, crit- 
ical, and supercritical linear-fractional BGW-processes with countably many 
types in the positively recurrent case. Finally, in Section 11 we summarize 



the implications of Section 10 results to the linear-fractional CMJ-processes. 

All our results for the linear-fractional BGW-processes with parameters 
(H, g, m) also apply to the case of finitely many, say p, types after putting 
zf"* = 0, = for 2 > p + 1, and hij = for i G [l,p], j > p + 1. The 
transient and null recurrent cases will be addressed in a separate paper. 



2. Linear-fractional distributions 

We are using the following notation for the vectors x = (xi,X2, . . .) of 
infinite dimension: 

= (0, 0, . . .), 1 = (1, 1, . . .), e, = (l|,=i}, 1|,=2}, . . .), 

the transpose of the vector x will be denoted as x*, and I will stand for 
the unit matrix (l{i=j})j>ij>i- Recall that for a non-negative random vector 
Z = (Zi, ^2, . . .) with integer-valued components its probability generating 
function is defined by E(s^) := ^.{sf^s^^ ...). 
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Definition 1. Let {ho, hi, h2, . . .) be a probability distribution on {0, 1,2,.. .} 
and {gi, g2, ■ ■ ■) be a probability distribution on {1, 2, . . .}. For a given positive 
constant m the function 

0(s) = ho+ , ^-^^t , (3) 

generates a probability distribution to be called a linear-fractional distribution 
LF{h,g,m), where h= (/ii, /i2, • • g= (fi-i, 5-2, • • •)■ 

Slightly modifying Theorem 1 from |12j (devoted to the finite-dimensional 
case) one can demonstrate that for a linear-fractional function 

J , , ao + aiSi + 02^2 + . . . 



bo + biSi + 62S2 + • • • 
to generate a proper probability distribution it is necessary that this function 
can be written in the form (|3|. It means that Definition [T] covers all possible 
linear-fractional probability generating functions. 

Next we show that a multivariate linear-fractional distribution can be 
viewed as a multivariate geometric distribution modified at zero. Namely, we 
claim that a random vector Z with generating function ([s]) can be represented 
as 

Z = X+(Yi + ...+Y^)-l|x^o} (4) 

in terms of mutually independent random entities (X, A^, Yi, Y2, . . .). Here 
vectors X and Y^ have multivariate Bernoulli distributions 

P(X = ei) = hi, i> 1, P(X = 0) = ho, 
P(Y, = e,) = gi, t> 1, J > 1, 

and is a geometric random variable with distribution 

P(iV = k)= m'=(l + m)-''-\ A; = 0, 1, ... . 



Since 



E(s^-l{x^o}) = $^/i 

i=l 

00 



i=l 

1 

1 + m — ms^ 
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the linear-fractional formula ([s]) is equivalent to Q due to 
E(s^) = ho + E(s^ ■ l|x^o})E((sYO'^)- 

Representation Q suggests a simple two-dice-one-coin experiment gen- 
erating a random vector Z with a linear-fractional distribution LF(h, g,m). 
Firstly, an {ho, hi, . . .)-die is thrown once producing an outcome i E {0,1, . . .}. 
If i = 0, the experiment is stopped and the vector Z is assigned value 0. Oth- 
erwise, we put X = and continue the experiment. On the second stage 
we repeatedly throw an m-coin which shows heads with probability and 
tails with probability Let be the number of tails before we see heads 
for the first time. Finally, a {gi,g2, . . .)-die is thrown N times to find the 
values of Y^. It remains to use Q to compute the corresponding value of Z. 

Definition 2. Let H = {hij)°°^^ be a sub-stochastic matrix with rows hj = 
{hii, hi2, ■ ■ ■) having non-negative elements such that hiQ := 1 — hn — hi2 — . . . 
take values in [0, 1]. Let {gi, g2, . . .) be a probability distribution on {1, 2, . . .}, 
and m be a positive constant. A branching process will be called linear- 
fractional with parameters (H,g,m), if its particles of type i reproduce ac- 
cording to the LF{hi, g, m) distribution, i = 1,2, . . .. 

Notice the strong limitation on the reproduction law requiring parameters 
(g, m) to be ignorant of mother's type. This is needed for iterations of the 
corresponding generating functions 

0;")(s)=E(s^'"'|zW = e,), z = l,2,... 

to be also linear- fractional. It is easy to see that if the denominators in 

0i(s) = hiQ + — , 

1 + m-ml^-^^gjSj 

were different for different i, then the iterations of these generating functions 
0j(0i(s), 02(s) . . .) would lose the linear-fractional property. 

3. Propagation of the linear-fractional property 

The matrix ^ of the mean offspring numbers for a linear-fractional 
branching process with parameters (H, g, m) is computed as 

M = H + mHl*g. (5) 
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After multiplying this by 1* we obtain 

Ml* = (l + m)Hl\ 
which leads to a useful reverse expression of H in terms of M 

TTl 

H = M Ml*g. (6) 

1 + m 

Theorem 3. Consider a linear-fractional branching process with parameters 
(H, g, m) starting from a type i particle. Its n-th generation size vector Z^"'^ 
has a linear- fractional distribution LF{h.["'\g^^\m^'"'^) whose parameters sat- 
isfy 



n-l 
1=0 

m(")g(") =mg(I + M + --- + M"~^), (8) 

(n) 

HW = m" -M"l*g("), (9) 

1 + m'^^' 

where Mi are given by ([T]) and H*^") is the matrix with the rows (h-"^)^]^. In 
particular, 

P(Z(") ^ 0) = (1 + m("))"iM"l\ (10) 
where P(Z(") ^ 0) is a column vector with elements P(Z'^"'' ^ 0\Zi^^^ = ej). 



Proof. Our proof of the fact that the ra-th generation sizes have a linear- 
fractional joint distribution is based on the contour process method and is 
postponed until Section |5| 

Turning to the proof of relations ([t]), ([s]), (|9|, and (10), observe that 
relation ([9]) is a straightforward counterpart of ^6]). The left hand side in 
equals H^^^l* so that (|To|) follows from (g. The only remaining relations 
to prove are ([T]) and ([s]). We do this using the spinal representation of the 
BGW-tree illustrated in Figure [Tj 

Both the spinal representation and the contour process methods (the 
latter discussed in Section |5| rely on a planar genealogical tree connecting the 
particles of the branching processes appeared up to the time of observation. 
For the current setting of linear-fractional branching processes with countably 
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Figure 1: A spinal representation of a BGW-tree reaching the observation level, cf [?]• 



many types it is important to use a particular plane version of the genealogical 
tree: given a group of siblings stemming from the same particle the leftmost 
branch should connect the mother to its first daughter (that one whose type 
may depend on mother's type). 

It is convenient to start with the proof of equality ([T]) even though it is an 
immediate consequence of ([s]) and ([T]). Suppose that at the observation time 
n population is not empty. The corresponding spine of the planar BGW- 
tree is the leftmost lineage of particles reaching the level n. The number 
of all branches present at level n except the spinal one has the mean m^'^\ 
Equality ([T]) simply says that this mean is the sum of contributions from all 
the lineages stemming from the spine to the right of it, see Figure [TJ In the 
linear-fractional case at each level / G [0, n — 1] there is a geometric with mean 
m number of branches growing off the spine to the right of it. Every one of 
such daughter branching processes, in accordance with ([T]), gives on average 
Mn-i-i descendants in generation n. Summing the products mMn-i-i over 
all spinal particles at levels / = 0, . . . , — 1 we get ([T]). 

Equality ^ is just a detailed version of ([T]) taking into account the num- 
bers of particles of various types existing at time n. 

□ 

Corollary 4. Conditionally on non- extinction T}-^^ has a multivariate shifted 
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geometric distribution 



E[sZ^"'|ZW^O,zW = e, 



y-oo ,(n) 



1 + mW -mW 



where W^ = f^^/{l-h^). 



Proof. The asserted formula follows from 
E[s^'"'l{z(")^o}|Z(°) = e,] = E[s^<"'|ZW = e,] - P[ZW = 0|ZW = e,] (11) 



after dividing the left and right sides by 



□ 



Example 5. As a simple illustration take a single type BGW-process with 
a shifted geometric offspring distribution 

= m'^-l(l + m)-^ A; = 1,2,.... (12) 

Clearly, this is a supercritical process with the mean offspring number 1 + m. 
The following algorithm introduces a labeling system that will make this 
branching process look like a branching process with infinitely many types: 

• the initial particle born at time zero gets label 1, 

• for any particle of type i its first offspring is labelled i + 1 while the 
remaining offspring are assigned label 1. 

The labeled process at time n will contain particles of types 1 through n + 1 so 
that the vector Z*^") of sub-population sizes will have components Z^""^ > 
for i = l,...,n, ^^+1 = 1, and Z^"'^ = for all i > n + 2. This is a 
linear-fractional branching process with parameters (H,ei,m), where H = 
{l{j=i+i})f^j=i, see Definition [2] 



9 



In this case it is easy to compute the mean numbers of particles at any 
time 



/ m(l + m)"~^ 
m(l + m)"~^ 
m(l + m)""^ 
m(l + m)""^ 

V 



m(l + m) m 1 

m(l + m) m 1 

m(l + m) m 1 

m(l + m) m 1 



\ 



and according to ([T]) we get 

i-i 



(n) 




m] 



n-1 

m^(l + 

j=0 



m] 



:i + m)"-l, 



which is verified by the fact that the mean of the n-th generation size in the 
original branching process is (1 + m)" = 1 + m^"^. This and (IS]) yield 



An) 



m 



n-1 



m 



(n) 



1 + m(l + my 



(1 + m)" - 1 



for i = 1, . . . ,n and g^"'^ = for i > n + 1. Finally, using ^ we obtain 

Combining these results we conclude that the components of the vec- 
tor , 2^^"'') are independent geometric random variables with means 



(m(l + m 



,n-l 



m(l + m 



in-2 



4. Individuals evolving in the type space 

A BGW-process is a population model of reproducing particles with non- 
overlapping generations. Here we give another interpretation of this popula- 
tion model in terms of individuals with overlapping generations. This view 
becomes very useful in the multi-type linear-fractional case. 

We start by explaining what we view as the initial individual for a given 
realization of the underlying BGW-tree. The initial individual is represented 
by the starting particle, its first daughter if any, the first daughter's first 
daughter if any, and so on. Such an individual has a random life length 
and it dies when a particle in the chain of the first-born descendants fails 
to produce offspring. In the multi-type setting the individual evolves in the 
type space as the particles forming the individual change their types in time. 
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Each moment of its life the initial individual is able to produce new 
individuals also evolving in the type space as sequences of the first-born 
particles. Turn for visual help to Figure [2jC giving the individual based 
picture of the same genealogical tree as in Figure [2j A. Each vertical arrowed 
branch in Figure [2] C represents an individual giving birth possibly multiple 
times during its life. In particular, the initial individual lives two units of 
time producing two daughters: one of them lives two units of time and the 
other only one. Notice that the first granddaughter of the initial individual 
produces two daughters at different ages. 

The evolution of an individual can be modeled by a Markov chain whose 
state space {0, 1, 2, . . .} is our type space {1,2,...} enhanced with a graveyard 
state {0}. The corresponding transition probabilities are given by the sub- 
stochastic matrix H together with hoj = for j > 1 and /iqo = 1- In terms 
this Markov chain the random lifespan L of a new-born individual is the time 
to visit the graveyard state starting form a state chosen at random among 
{1,2,...} according to distribution g. Note that the resulting life length 
distribution belongs to the broad class of discrete phase-type distributions, 
see for example [Ij. 

The above given description of the life length distribution implies that 
the tail probabilities for L are computed as 

dn := F{L >n) = gH"l*, n>0, (13) 

and the distribution of the random type of an individual at age n is given 
by the vector gH"/(i„. The latter observation brings the following natural 
condition which excludes the phantom (never observed) types. Let S be the 
set of types ever observed in the linear-fractional branching process, that is 
j G S, if for some n > 1 the j-th element of gH" is positive. Let p = \S\ be 
the number of elements in 5*. We can always relabel the non-phantom types 
in such a way that 

S = {l,2,...,p}zfl<p<oo, orS = {l,2,...} z/p = oo. (14) 



The generating function of the tail probabilities ( 13 ) 



f{s) = J2dns'' (15) 



n=l 



plays an important role in our analysis. Under condition (14) the function 



f{s) given by (15) is strictly monotonely growing from zero to infinity as s 



goes from zero to infinity. 
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Definition 6. Let 

Rf = inf{s > : f{s) = oo}. 

Define a positive finite constant R as R = Rf if f{Rf) < 1/m, or else, if 
f{Rf) ^ 1/"^, as the unique positive solution of the equation 

mf{R) = 1. (16) 

For such defined value of R set 



(17) 



n=l 



Lemma 7. Assume that condition (14) is valid. The following chain of im- 
plications is true 



{R<1} C{R< Rf} C {/3 < oo} n (16). 



Proof. Since dn < 1, we have /(s) < oo for all s G [0,1) implying Rf > 
1 and the first implication. The second implication is a straightforward 
corollary of Definition [6j 

□ 



To continue with important properties of the function ( 15 ) observe that 

/(■s) = ^~i-s ^ ~ ■'- ™plyiiig that the life length has always a positive mean 
(finite or infinite) 

A:=E(L) = 1 + /(1). (18) 

At every age during its life (but not at the moment of death) each individual 
produces a geometric number of offspring with mean m. Thus the mean total 
offspring number per individual is 

/i = m(A-l) =m/(l). (19) 



Observe also that given (16) 



f{Rs) = J2 dns"", L = mR^dn, (20) 



n=l 
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is the generating function for a probability distribution for a random variable 
(called the regeneration age of the immortal individual in [TT]) with mean /3. 
Notice that in the critical case dn = P(-Z^ > n)/{E,{L) — 1). 

Example 8. Observe that for any sequence 1 = > di > d2 > ■ ■ ■ with 
(i„ > for all n > 1 one can find a suitable pair (H,g) such that (13) holds. 
Indeed, choosing 



H 



/ 





V ••• 



rfi 

di/d2 

da/rfs 

d^/di 



\ 



J 



g = ei, 



we obtain gH"l* = dn- This observation allows us to illustrate the possibility 
of /3 = oo by letting d^ ~ cn~^r~'^ as n — >■ oo for some constants r > 1, c > 0, 
and A; > 0. Clearly, in this case Rf = r and 



if < /c < 1, then (16) holds with R <r and (3 < oo. 



if A; > 1 and /(r) < 1/m, then R = r and (16) does not hold. 



if A; > 1 and /(r) = 1/m, then (16) holds with R 



if A; > 1 and /(r) > 1/m, then (16) holds with R < r, 



if 1 < A; < 2, then /3 = oo, while for A; > 2 we always get /3 < oo. 



5. Jumping contour processes in discrete time 

Recall that for any planar tree, one can define a contour profile of the tree 
by the depth-first search procedure. Figure [2] illustrates the basic definition 
of the contour process for a finite tree supplied with a path around the tree. 
The contour process is simply the seesaw line graph (panel B) representing 
the height of the location of an imaginary car driving with a constant speed 
along the path outlined on the panel A (the imaginary car stops when it 
reaches the level —1). If the realization of the genealogical tree is infinite, 
one still can work with the contour processes after cutting off the branches 
above the level n corresponding to the observation time, as demonstrated in 
Figure [2] It is easy to reconstruct the tree from the contour process. For 
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Figure 2: A: a BGW-tree (thick line) stopped at level n = 5 and supplied with a contour 
(dashed) line. B: the corresponding unfolded contour process. C: the CMJ-view of the 
same tree when depicted in terms of individuals forming a branching process with overlap- 
ping generations. The vertices marked by arrows represent individuals which are dead at 
that time. The stopped tree gives no information about the fate of the three tip vertices. 
D: the modified contour process of a constant speed descent with iid upward jumps. To 
match the figure on panel B the jumping contour process starts at the level -1. 

example, the number of particles alive at time n, if any, is 1 plus the number 
of excursions of the contour process starting at level n downwards and coming 
back to the level n escaping absorption at level -1. In particular, in Figure [2] 
there are two such excursions for n = 5 with the number of branches at this 
level being 3. 

The contour process approach has proven to be very useful in the theory 
of branching processes (see for example [,8] and references therein). In the 
single type linear-fractional case the contour process has a very simple struc- 
ture of an alternating random walk. Alternating upwards and downwards 
stretches have independent lengths following two shifted geometric laws (one 
for upward and the other for downward moves). 

In the multi-type linear-fractional setting, one can ensure a Markov prop- 
erty of the contour process using additional labeling so that the states of 
the contour process will be given by pairs (Z,i) with I = —1,0, 1,2, .. . and 
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i = 0,1,2,.... The current state with i > I tells three things about 
the contour process: it is visiting level /, the last move was upwards, and 
the corresponding particle in the branching process is of type i. The states 
(/, 0) are attained by the contour process after the downward steps. This 
convention implies in particular, that at the absorption level I = —1 the only 
possible label is z = 0. Clearly, the labeled contour process is a Markov chain 
with transition probabilities formulae 

P{(/,z) -> (/ + l,j)} = h^,, ^ (/ - 1,0)} = /i,o, 

TTl 1 

P{(/, 0) (/ + 1, j)} = y—9j, 0) ^ (/ - 1, 0)} = -— , 

1 + m 1 + m 

P{(-1,0)^(-1,0)} = 1, 

vahd for all z > 1, j > 1, / > 0. 

The following alternative way of introducing Markovian structure in the 
contour process of a linear-fractional multi-type branching process does not 
require additional labeling. What we call here the jumping contour process (cf 
[I7j) has a trajectory of a constant speed descent with independent upward 
jumps each distributed as the individual life length L. From any given current 
level / the jumping contour process moves one level down to / — 1 and either 
settles there with probability or, with probability jr^, it instantaneously 
jumps L levels up coming to the level l — l + L. Figure |2jD clearly illustrates 
the last construction. 

Proof of the first statement of Theorem [3l Consider now a linear- 
fractional multi-type BGW-tree stopped at the level n so that the set of its 
labeled tips allows to compute the current generation sizes Zi^"'\ Using the 
contour process idea it is straightforward to see that the distribution of the 
vector Z*^"^ must be linear-fractional. 
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We have to verify that if there is at least one branch reaching the level 
n, then the number of tips and their types at the level n to the right of 
the leftmost tip has a multivariate geometric distribution. In terms of the 
jumping contour process the consecutive tips to the right of the leftmost 
tip are connected by a geometric number of independent excursions starting 
from the observation level n and returning from below back to the level n. 
The parameter of the corresponding geometric distribution is the probability 
that the jumping contour process starting from level n will be absorbed at 
the level -1 without visiting the level n once again. This finishes the proof of 
the first statement of Theorem [3l 



6. Classification of branching processes with countably many types 

Branching processes are classified according to the asymptotic properties 



as n oo. The assumed independence of particles implies a recursion 
M(") = MM("-i), where M = M^^). It follows that M^") = M". Given that 
all powers M" are element-wise finite (which is always true in the linear- 
fractional case) the asymptotic behavior of these powers is described by the 
Perron- Frobenius theory for countable matrices (see Chapter 6 in [22]). 

Next we remind some crucial conclusions from this theory holding for an 
irreducible and aperiodic countable matrix M. Recall that a non-negative 
matrix M is called irreducible, if for any pair of indices (z, j) there is a natural 
number n such that m[^^ > 0. The period of an index i in an irreducible 
matrix M is defined as the greatest common divisor of all natural numbers 
n such that m,-"'* > 0. In the irreducible case all such indices have the same 
period which is called the period of M. When this period equals one the 
matrix M is called aperiodic. 

Due to Theorem 6.1 from [22] all elements of the matrix power series 




□ 





oo 




(21) 



n=0 
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have a common convergence radius < i? < oo, called the convergence 
parameter of the matrix M. Furthermore, one of the two alternatives holds: 

oo 

/^-transient case: m[^^R^ < oo, z = 1, 2, . . . , (22) 

n=0 

oo 

i?-recurrent case: mff^ B!^ = oo, i = 1, 2, . . . . (23) 

n=0 

According to [22] (Theorem 6.2 and a remark afterwards) in the i?-recurrent 
case there exist unique up to constant multipliers positive vectors u and v 
such that 

i?Mu* = u*, i?vM = V. 

A renormalization Rvjrriji/vi transforms the matrix M into a stochastic ma- 
trix. 

The i?-recurrent case is further divided in two sub-cases: i?-null, when 
vu* = oo, and i?-positive with vu* < oo. In the i?-null case (and clearly also 
in the /^-transient case) 

i?"m;;)^0, t,j = l,2,... 

In the -R-positive case (Theorem 6.5 from [22]) one can scale the eigenvectors 
so that vu* = 1 and obtain 

R^m^^ ^u,v„ z,j = l,2,.... (24) 

These results suggest a double classification of the branching processes 
with countably many types having a mean matrix M. The usual classification 
depends on the value of the Perron- Frobenius eigenvalue p = 1/R. Given 
p < 1, p = 1, or p > 1 the branching process is called subcritical, critical, 
or supercritical. An additional classification is needed to distinguish among 
different asymptotic patterns due to recurrence property of the branching 
process with respect to the infinite type space. 

Definition 9. A branching process with countably many types will be called 
subcritical (critical, supercritical) and transient {recurrent, null-recurrent, 
positively recurrent} in the type space, if its matrix of the mean offspring 
numbers M has a convergence radius R > 1 (R = 1, R < 1) being R- 
transient {R-recurrent, R-null recurrent, R-positively recurrent}. 
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There are several published results for the BGW-processes with countably 
many types (see for example [2], [15]). One of them is Theorem 1 in [18] 
dealing with the i?-positively recurrent supercritical {R < 1) case. It states 
that if 

oo 

^t;,E((Z«u')2|z(°) = e,) < oo, 

i=l 

then for any w such that w < cu for some positive constat c, the convergence 
i?"Z„w* — Fvw* holds in mean square, where Y > has a finite second 
moment. This statement is cited here just to illustrate the need for finding 
examples of branching processes, where conditions like -R-positive recurrence 
could be verified and the values of {R, u, v) be computed in terms of the 
basic model parameters. The linear-fractional framework allows for such 
amendments. 



7. Main result 

Consider a linear-fractional branching process with parameters (H, g, m) 
whose matrix of the mean offspring numbers M is given by ([s]). Clearly, if 
M is irreducible, then starting at any type it is possible to observe any type 



in S at some later time. Thus irreducibility of M implies (14). The opposite 
is not true, if there exist so-called final types, that never produce offspring. 
To exclude the final types we will assume that 

there are no zero rows in H. (25) 

This is of course also a necessary condition for irreducibility of M. 

Before stating our main result we remind that according to Lemma [7] 
parameter (3 defined by ([17]) can be infinite only when i? > 1, see Definition 
|6} Due to the same Lemma [T] condition f{R) < 1/m is only possible again 
when R>\. 

Theorem 10. The matrix M given by ([s]) is irreducible if and only if (14) 
and (25) hold. IflSA is irreducible and aperiodic, the following five statements 



are valid: 



(i) the convergence parameter R o/M is computed from the gen- 
erating function (15) as described by Definition^ 



(ii) if f{R) < 1/m, then M is R-transient, 
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(iii) if f{R) = 1/m, then M is R-recurrent, 

(iv) if f{R) = 1/m and (3 = oo, then M is R-null recurrent, so 
that each element of R^M."^ converges to zero as n ^ oo, 

(v) if f{R) = 1/m and (3 < oo, then M is R-positively recurrent 
so that 

R^'M'' uV, n oo, (26) 

where element-wise positive and finite vectors u and v are 
given by 



u* = (l + m)/3-' 

k=l 

oo 



m 



1 + m 



(27) 
(28) 



satisfying vu = vl = 1 and gu 



fc=0 

t 1 r<-.it l+m 



mh 



Example 11. Consider a triplet (H, g,m) such that gH = rg for some 
positive constant r. Since H is sub-stochastic, we have necessarily r < 1: 

r = rgl* = gHl' < gl* = 1. 

It follows from (|5| that g is a left eigenvector of the matrix M as well: 
gM = r(l + m)g. Using (16), (17), and (28) we obtain 



p = (1 + m)r, (3 



l+m 



m 



Using (13) we obtain P(L > n) = r". In this case the branching process is 
always i?-positively recurrent. 

Example 12. Assume now that HI* = rl*, where necessarily r < 1. It 
follows from (jsj that Ml* = r(l + m)l*. In this case 



p = (1 + m)r, (3 



m 



u = 1. 



m 



It follows from (13) that even in this case P(L > n) = r". The case with 
r = 1, giving the infinite life length as in Example 5 from Section [3| is ob- 
tained whenever the matrix H is stochastic. 
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Remark 13. In the linear-fractional setting we have two alternative forms for 
the above given criterium for criticality recognizing three classes according to 
p < 1, p = 1, and p > 1. In terms of the mean number of offspring p produced 



by a single individual, see (19), the linear-fractional BGW-processes with 
countably many types can be called subcritical, critical, or supercritical if 
p<l,p = lorp>l respectively. Observe that in the i?-transient case 
with p = 1 we get p < 1, since /(I) = f{R) < l/m. 

On the other hand, super-criticality can be identified via a positive drift 
for the contour process. The drift A — 1 — of the jumping contour process 
is computed as the difference between the mean jump size L and the length 
of a downward stretch 1 + m~^. Clearly, the inequality A — 1 — m^^ > is 
equivalent to p > 1. 

8. Renewal theory argument 

This section contains two lemmata which will be used in the next section 



proving Theorem 10 The following well-known renewal theorem taken from 



Chapter XIII. 4 in [5j will be used by us several times. 

Lemma 14. Let A{s) = Yl'^=o^nS^ be a probability generating function and 
^{^) — Xl^o^"'^" ^'^ ^ generating function for a non-negative sequence so 
that A{1) = 1 while B{1) G (0,oo). Then the non-negative sequence defined 

by Z]^r=o^"'^" = i^W) ^^^^ ^^"^ ^" ~^ as n ^ oo. 



The next key lemma deals with the matrix power series (21). 



Lemma 15. Let R be given by Definition^ For any s G [0,R) the matrix 
power series (21) satisfies 



777 

M{s) = H{s) + Y3^^(H(.) - I)1%H(.), (29) 

where H(s) = X]^o ■^"H'^ is element-wise finite. Moreover, 

gM(i?)l*<oo, z//(i?)<l/m, (30) 

gM(i?)l* = oo, z//(i?) = l/m, (31) 

gM(s)l* = oo, for s> R. (32) 
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Proof. Due to Q we have M"+i = (H + mHG)M", where G = l*g. 
Using induction we obtain 



M" = H" + m ^ H*GM"-\ 

i=l 

This yields a renewal equation for the sequence M„ = gM"!* 

n 

Mn = dn + m'Y diMn-i 

1=1 

which translates in terms of generating functions into 

oo oo 

Mns^ = 1 + f{s) + mf{s) Y Mns""- 



(33) 



n=0 



n=0 



From here it follows that for s G [0, R) the generating function 



gM(s)l* 



1 - m/(s) 



(34) 



takes finite values, and that relations (30), (31), and (32) hold. 

Let s G [0,-R). Relation (34) entails that M(s) can not be element-wise 
infinite and therefore, the radius of convergence for M is larger or equal to 
the value R given by Definition |6} Therefore, relation (33) implies 

M(s) = H(s) + msH(s)GM(s) 

= H(s) + msH(s)GH(s) + . . . + (msH(s)G)"M(s) 

with all terms being element-wise finite. Observe next that for all > 1 

/ oo N 

(sH(s)G)"= K^sWl*g 



oo 



.1=1 



Y s'H'l' Y ) • • • I XI ^'gH^l* ) g 

Ki=l J \i=l 

OO 

r-i(.)5^.wi*g. 



i=l 
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Thus due to mf{s) < 1, the term 

oo 

3H(s)G)"M(s) = mT-\s) ^ s^ffl*gM(s) 



[msi 



vanishes as n — > oo and the previous two relations yield (29) 

oo 

M{s) = J](sH(.)A)"H(.) 



n=0 

oo oo 

His) + J2 ^"T-'is) J2 sWgiiis) 

n=l 1=1 

(U(s)-I)l's-U(s). 



n{s) + - 



m 



mf{s) 



□ 



9. Proof of Theorem 



We show first that conditions (14) and (25) are sufficient for irreducibility. 
Namely, we will demonstrate that for any j there is a n = rij such that 



m^"^ > for all i. To see this observe that (33) gives 



M" > mH^l*gH"- 



for any k G [l,n]. Due to (25) for any k the i-th component of the vector 
H'^I* is positive. On the other hand, under (14) for the given j there is an / 



such that the j-th component of the vector gH is positive. 

Assume from now on that M is irreducible and aperiodic. Statement (i) 
of Theorem 10 follows from Lemma 15 as M(s) is element-wise finite for 



s G [0, R) and must be infinite for s > i? due to (32). 

Statement (ii) is an immediate consequence of (30): given f{R) < l/m, 



relation (23) is not possible, and we are in the transient case (22). 



To prove (iii) we assume that f{R) = l/m and verify that M(i?) is 
element-wise infinite again using Lemma 15 Indeed, since H(s) > I, each 



element of the vector gH(s) is larger or equal than 1, implying according to 
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(29) that for any pair of indices and any s G [0, R) 



— mf[s 

Ci{s) 
1 - m/(s) 



where Cj(s) is strictly positive due to (25). Letting here s — )■ i? we see that 
every element of M(i?) is infinite as soon as f{R) = 1/m. 

Suppose ( 16 ) holds. Assertion (iv) is easily obtained using Lemma 14 It 
follows from^34|) that i?"gM"l* given /3 = oo. Thus (|24|) can not be 
true in this case implying null- recurrence. 

It remains to prove (v). Let /3 < oo. We show next that H(i?) is element- 
wise finite. To see this notice that gH(i?)l* = 1/m according to (16). On 
the other hand, under condition (14) for any index i there is k = ki and a 
positive Ci such that gH'^ > implying 



ra=0 



n=0 



and further 

eiH(i?)l* < c-^R-^m-^ < oo. 

Consider vectors u and v which, thanks to the just proved finiteness of 
H(i?), are well-defined by (27) and (28). The claimed equality vu* = 1 
follows from 

oo 

mgH{R){H{R) - 1)1* = m E ni?"gH"l* = /3, 

n=l 

which is a consequence of 

oo oo 

H(s)H(s) - H(s) = J2 H's^ J2 ^^^'^ 

1=1 k=0 

oo oo oo n oo 



1=1 n=i 



n=l i=l 



n=l 
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Finally, to prove ( 26 ) define a sequence of matrices B„ by 

oo 

m 



n=0 



■mf{Rs) 



[H(fe) -I)rgH(i?s) 



According to Lemma 14 we have an element-wise convergence B„ — )■ u*v as 



n — )■ oo. It remains to see that due to (29) 



where each element of K^li"' converges to zero, since H(i?) is element- wise 
finite. 



10. i?-positively recurrent case 



Propositions 16, ^71 and 18 below deal with linear-fractional branching 



processes given that M irreducible and aperiodic and that (16) holds with 



(3 < oo. Due to Theorem[TO]we are in the i?-positively recurrent case, so that 
M" ~ p"u*v, where p = R~^. Notice that the left eigenvector v describes 
the stable type distribution: ejM" ~ Uip"'\, and the right eigenvector u 
compares productivity of different types: M"l* ~ p"u* (so that Ui can be 
interpreted as the "reproductive value" of type i). 

These propositions should be compared with the known versions of such 
statements in the finite-dimensional case given in [12] and [20]. We remind 
about formula ( [sj) a vailable for M, formula (17) for (3, (27) for u, (28) for v, 
(pisl) for A, and PI for p. 



Proposition 16. In the subcritical positively recurrent case when p < 1, or 
equivalently p < I, 

P[Z(") 7^ 0] ~ p"(l + m)-\l - p)u\ (35) 
and independently from the initial type i we have a convergence in distribution 

[zW|z("Vo,zW = e,] Ay 

to a random vector Y with a linear-fractional distribution LF{h,g,m), where 

h = (1 + m)(l - - mg(I - My\ 

g = \'\l-p)g(l-M)-\ 
m = mA(l — p)^^, 



with hi* = 1. In particular, Yl* has distribution (12). 
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Proof. From ^ and ([34]) we obtain 



1 + 

l-m/(l) 



which together with (10) imphes (35) . The statement on the convergence of 
the conditional distribution of Z'^'^foUows from Corollary |4j 



1 + m -m}_^j^^gjSj 

since m'^^^g*^") — )■ mg(I — M)^-*^ = mg and 

H(") ~ p"(uV - (1 - + m)-^u'mg) = p"(l - + m)-^u'h. 



□ 



Proposition 17. /n i/ie critical positively recurrent case when p = 1 we have 

P(Z('^VO)~ n~\l + m)-^/3u\ 

and 

[n-iz(")|Z("V0,Z(°) = e,] 4 Xv, 

where X is exponentially distributed with mean (1 + m)P^^. 
If furthermore a vector w is such that vw* = 0, then 



[^-i/22(")w*|z(") ^ 0, Z(°) = e^] A rv^(l +m)(2/3)-iw2vt, (36) 

where = (w^, w|, . . .) and Y has a Laplace (double exponential) distribu- 
tion with parameter 1. 



Proof. Lemma 14 and relations ([t]), ([34]) imply that in the critical case 

m^") ~ n(l + m)/3-\ 
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Thus the stated asymptotics for the survival probabihty follows from (10). 
The second assertion readily follows from Corollary |4] as 



m 



and the multivariate shifted geometric distribution can be approximated by 
a one-dimensional exponential distribution along the vector v. 



Convergence (36) is verified in terms of the moment generating functions 

^ / (./^)z(")w*l2(n) / 2(0) _ \ _ E(sg'"'l|z(n)-,o||ZW =e, ^ 
V |Z- ^ u, - - ^ ^ 



[ZW ^ 0|Z(o) = e, 

where s„ = (e^^^^-'"'^ , e'-^''^^"'^ , . . .). Rewrite the numerator using (11) and 



observe that it is equal to (see Theorem [3]) 



In view of 

oo 

^ hf^e^'l^^^ ~ 1 - = P(Z(") ^ 0|Z(°) = e^), 

oo 

we conclude 

r 1^ ^"'^ V ^ 1-^2(1 + ^)(2/3)-iw2vf 

Proposition 18. /n the supercritical positively recurrent case when p > I 
P(Z("V0)^ {p-l){l + m)-^/3u\ 

and 

[p-"Z(")|Z(") 0, Z(°) = Bi] A Xv, 



□ 
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where X is exponentially distributed with mean (1 + m)(p — 1) ^ . 
If furthermore a vector w is such that vw* = 0, then 

[p-"/22H^t|2W ^ o,Z(°) = e,] A Y^/{1 + m){p - l)-l(2/3)-lw2v^ 

where = [wf, W2, ■ ■ ■) and Y has a Laplace distribution with parameter 1. 

Proof. Rewrite ^ as 

n-l 
1=0 

to obtain the following consequence of p4j) 



' il-mf{sR)){l-Rsy 



Thus Lemma [m entails 



m 



(n) 



p-{i + m)r\p-i)-'. 



This together with (26) and (10) gives the stated formula for the survival 
probability. The remaining two assertions on weak convergence are proved 
in a similar way as in the critical case above. 

□ 

Remark 19. To illuminate the above asymptotic results in the critical and 
supercritical cases leading to the Laplace distribution in the limit, we sketch 



a probabilistic argument for (36) in the two-type case when Z(") = (X„,iV„- 
Xn). Here conditionally on the total population size Nn the subpopulation 
sizes are asymptotically binomially distributed 

X„ ^ Bin(Ar„, vi), Nn-Xn^ Bin(Ar,, V2), 

and given WiVi + W2V2 = 0, we are interested in the asymptotic distribution 
of the linear combination 

Yn = WiXn + W2{Nn - X„) = {wi - W2){Xn - NnVi). 



27 



With a large the distribution of Yn/ y/n is asymptotically normal with 
zero mean and variance {Nn/n){wi — W2YV1V2 = A^„w^v*. On the other 
hand, conditionally on non-extinction, Nn/n is asymptotically exponentially 
distributed with parameter (1 + m)~^/3. The last two facts combined to- 
gether result in a limit Laplace distribution for Yn/ \fn obtained as a normal 
distribution with a random variance having an exponential distribution. 



11. The linear fractional CMJ-process 

In this final section we focus on the process formed by the total population 
sizes Z*^") = Z^")l*. As the particle types are neglected the process {-Z'^"^}J^q 
is not a Markov chain anymore, however, it gives an interesting example 
of a discrete time CMJ-process [H] which will be called a linear-fractional 
CMJ-process. Its individual life length, if unbounded, can have an arbitrary 
distribution, see Example 8, and the numbers of births given at different 
ages (except at the moments of birth and death) are independent geometric 
random numbers with the same mean m. 

For the linear-fractional CMJ-process to be fully described by the life 



length distribution (15) and parameter m one has to assume that the initial 
particle of the underlying linear-fractional branching process with parameters 
(H,g,m) has distribution g: P(Z^°) = e^) = gi, i > 1, so that the life 



length of the initial individual has the same distribution ( 15 ) as all new-born 
individuals appearing in the linear-fractional CMJ-process. 

In the framework of CMJ-processes (see |9|) the established classification 
of branching processes is expressed in terms of the so-called Malthusian pa- 
rameter a = —InR. It is assumed that R, a positive solution of equation 



(16) is well defined, a restriction which in the current framework of multi- 
type BGW-processes excludes from consideration the transient case. The 
usual criterium of criticality {a < 0, a = 0, a > 0} for the CMJ-processes, 
when a is defined, is of course equivalent to the classification {p < 1, p > 1, 
p > 1}. Another key parameter /3 is called the mean age at childbearing and 
plays an important role when the renewal theory methods are applied for the 
analysis of the CMJ-processes. 

Next we give a summary of the results for the linear-fractional CMJ- 
process which follow from the more general results obtained earlier in this 
paper. In the case when the Maltusian parameter a exists and /3 < 00, the 
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asymptotics of the survival probability 

oo 

Q„ = P(Z(") > 0) = 5^ (7.P[Z(") ^ 0|Z(°) = e,] 



i=l 



according to Propositions 16, 17, 18 takes a particularly transparent form: 

g„~e""(l-/i)/(m/3), if a <0, 
Qn ~ l/(nm), if a = 0, 
Qn (e" - l)/m, if a > 0. 



Moreover, in the subcritical case due to Proposition 16 there exist a discrete 
limit distribution 

P(Z(") = > 0) ^ m^-^(l + m)-^ = 1, 2, . . . . 

In the critical case Propositions [17] implies that for any positive x 



Finally, in the supercritical case Propositions 18 yields a weak convergence 
result 

P(Z(") > e""x|Z(") > 0) -> e-'^^, X > 0, 

where a = /3(e° — 1)/(1 + m). These explicit results nicely illuminate much 
more general limit theorems obtained for the CMJ-processes, see for example 

[9], dn], and m- 



Remark 20. There are many triplets (H, g, m) resulting in the same linear- 
fractional CMJ-process {Z^"'^}. For all of them the values of R and (3 will 
be the same as they are governed by the individual life distirbution. The 
difference between such related branching processes is in the labeling rules for 
particles which can be seen only in the expressions for the Perron-Frobenius 
eigenvectors. 
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